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1. Introduction 

In [1,2] I computed the operator determinants on lens space factors of the 
unit 3-sphere for the Laplacian conformal in four dimensions. This meant that 
the eigenvalues were squares of integers and the (^-function analysis was relatively 
straightforward and explicit. Nash and O'Connor, [3], evaluated the determinants 
for other operators in connection with analytic torsion. Factors of higher spheres 
have been considered by Bauer and Furutani [4]. Evaluations have also been con- 
ducted lately in string theory and quantum field theory, e.g. Castro, Lashkar and 
Maloney, [5], Gang, [6], Alday, Fluder and Sparks, [7], Radicevic, [8]. 

I here return to this topic still restricting to the 3-sphere, essentially because 
it is somewhat special and therefore more interesting (and easier). My attitude is 
mainly numerical and I wish to generalise the Laplacian to one that includes mass, 
in particular to the Laplacian conformal on the three-sphere. The analysis may be 
useful when considering the topology of the Universe, a topic under present scrutiny. 
I feel that a mostly computational approach has advantages such as speed and will 
prove a useful adjunct to the analytical route. 

Fixed-point-free factors of the three sphere, S 3 /r', divide into the homoge- 
neous variety in which the deck group, T', acts on the right (or left) the left (or 
right) action being trivial, and the inhomogeneous sort with a two-sided action. 
Left and right here refer to the SU(2) groups in the symmetry group isomorphism 
SO(4) ~ SU{2) x SU(2)/Z 2 acting on S 3 in its guise as SU(2), i.e. V = T' L x T' R . 

In this paper I mostly confine my attention to the homogeneous type which is 
somewhat easier to deal with. The classification of relevance here is given in Wolf, 
[9], Corollary 2.7.2 which says that T' R is either a cyclic group or a binary polyhedral 
group 2 . 

I showed in [1] that, because of the homomorphism SO(3)~ SU(2)/Z 2 , the 
binary degeneracies, which are the essential part of the spectrum, can be reduced 
to an analysis on the orbifolded two-sphere, S 2 /r, where V is the lift of V. Following 
Bethe, I refer to V as the double of V. For many purposes, one does not need the 
degeneracies themselves, their generating functions often being sufficient, and these 
are classically available, e.g. [11]. 

Lens spaces arise from the choice of a cyclic group, Z , for r^. In order to use 
the reduction mentioned in the last paragraph, the restriction that q should be even 

2 I consider the dihedron to be a polyhedron. Incidentally, my favourite potted summary of the 
classification is that by Ellis, [10]. 
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has to be made because the binary lift of Z ? is Z 2? . However it is possible to treat 
even and odd orders together, which I will do. 

2. Massive determinants on homogeneous lens spaces 

In order to be in accord with a previous notation, the scalar eigenvalues are 
taken to be I 2 — a 2 , Z = 1,2,..., with total degeneracies, D^V). I need the associated 
^-function which I write using a Bessel function relation, [12,13], as, cf [14], 



if/°°*EA(r')^(^)'" 1/2 /,- 1/2 (-) 



(1) 



-1/2 . 



where K r , is the cylinder kernel for the Laplacian conformal in four dimensions, 
and is, equivalently, the degeneracy generating function used in [1]. 

For the homogeneous case, the total degeneracy is a product of the left and 
right degeneracies, the former simply being a factor of /. (One should regard / as 
an SU(2) representation dimension, 2j + 1, which is just the range of the left index 
on a Wigner matrix.) Then, D^T') = Id^T') and we need the second factor. 

For quotient spaces, simple character theory allows a closed form for the gen- 
erating function of right degeneracies. In particular, for lens space degeneracies, 
d t {q), 

G(t,,) = f>(^ = j^J^-y ( 2 ) 

and, therefore, setting t = e _T , 

oo 

Kl' 2 (r) = Y,ldM)e- W 

l=i (3) 
d coth(gr/2) 

dr 2 sinh r 

The mathematical problem before us is the evaluation of the derivative at zero, 
Z'(0,q, a), of the (-function (1) with (3). Since the dimension of the manifold 
is odd, I can proceed exactly as in [13], [14], [15] and employ a complex contour 
method of continuing to s = 0. 
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Defining I(s, q, a) by Z(s, q, a) = I(s, q, a) /T(s), one has Z'(0, q, a) = 1(0, q, a) 
and a simple calculation which parallels that in [14] produces, 



1(0, q, a) 



o 



00 (q sinh r + cosh r sinh qr) cosh ar 

\JbJb JTVt; " ^ 2 " 



2rsinh r sinh qr/2 



(4) 



where r = x + iA with A < 2iv/q. For S 3 , q = 1 and 
7(0, 1, a) = 



o 



00 (1 + cosh r) cosh ar 
ax Ke : 2 



Re 



2rsinhr sinh r/2 
cothr/2 coshar 
2r sinh 2 r/2 ' 



which, as a check, agrees with an expression derived in [14]. 

In Fig.l, I plot logdet = —1(0, q, a) as a continuous function of q for a = and 
for a = 1/2, corresponding to conformal in four and in three dimensions respectively. 

Fig.1 . lens log det 
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Figure 2 gives, for three values of q, the variation of logdet against the "mass" , 
\i defined by /U 2 = 1/4 — a 2 , which is, for this purpose, a more convenient variable. 
(What should be called the mass is somewhat arbitrary.) 

The asymptotic variation for large mass was discussed, in a similar context, in 
[14] and is best expressed in terms of the short-time expansion coefficients of the 
heat-kernel for propagation by the Laplacian conformal in four dimensions because 
this expansion, in the present case, terminates with the first, 'volume' term and it 
easily follows that Z'(0,q,im) ~ nm 3 /3q with an exponentially small remainder. 
The beginnings of this can be seen in the figure. 
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Fig. 2. lens log det 
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3. The other quotients 

As a simple 'application' of these expressions, the logdets for the binary tetrahe- 
dral, (T'), octahedral, (O'), and icosahedral, (I'), factored spaces can be determined 
using the cyclic decomposition of the corresponding SO(3) groups, [11,16]. This was 
used in [1]. Any spectral quantity 3 S(V) on such an S 3 /r' can be expressed in terms 
of binary cyclic quantities, 

S{T') = \(Y, S ^2 q )-S{'L 2 )\ (5) 

\ q ' 

where the sum is over the (three) g-fold axes of rotation in the SO (3) subgroup, F, 
of which r' is the lift. The values of q are contained in the symbol of the polyhedral 
group, (I, m,n), i.e. I = 2, m = 3 and n = 3, 4, 5 for T, O, I, respectively. 

I give the numbers for the conformal-in-three and conformal-in-four dimen- 
sions Laplacians. The latter agree with those in [1] computed using a different 
algorithm for the even lens spaces. Thus, 

det (T') = 0.159259 , det (O') = 0.099650 , det (I') = 0.055743, 

3 I restrict the notion of spectral quantity by imposing linearity in the sense that, if the spectrum 
is composed of the union of two sets, its spectral quantity is the sum of the spectral quantities 
of the two parts. A typical example is the ^-function but not the determinant. 
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for conformal-in-three, and, 



det (T') = 0.202089 , det (O') = 0.128776 , det (I') = 0.073056, 
for conformal-in-four dimensions. These values are for real scalar fields. 

4. Twisted fields 

A simple, and basic, example is a complex scalar field or U(l) line bundle. In 
general terms, the quantum field theory is characterised by the homomorphisms of 
the deck group, r', which here is the fundamental group, into the bundle group, 
Q, i.e. by Hom(r',^). which can be realised as flat connections associated with 
'fluxes' of gauge fields through 'holes' in the configuration space, e.g. [17-20]. 

For the simplest situation, r' = 7L q and Q =U(1), as considered in [17], for 
example. The homomorphisms correspond to the q qth roots of unity and the 
elements of Hom(Z q , U(l)), i.e. the multiplying representations of Z q (the 'twists') 
are generated by u) r where u) q = 1 and < r < q. 

Instead of (2), character theory gives the degeneracy generating function 

G(t,q,r) = V d,(q,r) t l = - K^r 1 t, 

v ' y ' ' j-^ lKH " ' (1 -t 2 )(l -ti) 

where 8 + 1/2 = r/q, which is the flux. For (3) one has, 

oo 

K 1 q / 2 (r,r) = J2ld l (q,r)e- 1 - 

'=i (6) 
d cosh(qrS) 

dr 2sinhr sinh qr / '2 

The effect of the twisting on the logdet can now be determined numerically 
from the equivalent of (4), i.e. , (r = x + iA), 

l(0ara) = - f°° dx Re ( cos ^ 1 aT — cosh ^ Tj ) ^ ( 7 ) 

l{U,q,r,a) ^ ax tie ^ ^ dx 2smhr sinhqr/2 J ' UJ 

where I have not, this time, performed the differentiation. 

Figure 3 shows the (typical) variation of the effective action (minus half the 
logdet), for a complex field with the twist of the bundle plotted continuously. The 
mass parameter, a, has been set to 1/2 to give the Laplacian conformal in 3-d. 
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Fig. 3. twisted lens effective action 
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5. More complicated twistings 

More generally, if the field belongs to some representation of the internal sym- 
metry (bundle) group, Q, a twisting, i.e. an element p G Hom(r',(?) is determined 
by expressing this representation as a direct sum of irreps, A, of V . All that is 
needed, for a flat bundle, is that the dimensions match and that any other condi- 
tions, such as unimodularity, be met, [21]. This is an enumerative, rather than an 
algorithmic, prescription for cataloguing the gauge vacua. 

For a flat bundle spectral quantity, S(p), the direct sum becomes, by additivity, 
an algebraic sum of the components, 5(A), which act as spectral building blocks. 
These, in turn can be expressed linearly in terms of cyclic spectral quantities. This 
follows from the algebraic sufficiency of the cyclic subgroups of T', which is guar- 
anteed by Artin's theorem, and an inducing isospectral theorem of Ray and Singer, 
[22]. (See [23], Tsuchiya, [24].) Equation (5) is an example which I rewrite more 
explicitly as, 

5(1) = i (5(0; R) + 5(0; S) + 5(0; T) - 5(0; RST)) , (8) 

where 5(r;7) is the twisted spectral quantity for a cyclic quotient generated by 7. 
R, S and T are the generators for the binary polyhedral groups, e.g. [25], corre- 
sponding to q = 2, q = 3 and q = (3,4 or 5), respectively. 4 The left-hand side of 
(8) refers to the trivial irrep, 1, of T', i.e. the trivial bundle. 

Only R, S and T are really needed. The inclusion of the central element E = RST gives a 

2 

more symmetrical formulation, E = E = id. 
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As another example, derived in [23] , I give the first twisted relation, 

5(2.) = \ (5(1; R) + 5(1; S) + 5(1; T) - 5(1; RST)) , (9) 

where 2 S is the two-dimensional spinor representation of T', O' or I'. 
Numerical evaluation of (9) using (7) produces, 

det (T') = 0.652112 , det (O') = 0.439366 , det (I') = 0.260126, 

for conformal in three, and, 

det (T") = 0.663348 , det (O') = 0.454594 , det (I') = 0.272797, 

for conformal in four dimensions. 

Just to illustrate the sort of calculation that is possible, I consider a (scalar) 
field belonging to the fundamental rep of an internal U(3) on icosahedral space, 
S 3 /F. There are four elements of Hom(r,U(3)), viz 3, 3', 1 © 2 S , 1 © 2' 9 and the 
trivial one, 1 © 1 © 1. 5 Inducing representations shows that, [23], 

5(2 / J = 5(l;5)-^5(5;T)-5(l;T) 
5(3') = ^ 5(2; R) — 5(2; T) 
S{Z)= 1 -S{2-R)-S{A-T). 

Together with (8) and (9) these yield the following numbers for the 4d-conformal 
Laplacian, 

det(l© 1© 1) = 0.000391 
det (1 © 2J = 0.019929 
det (1 © 2' 8 ) = 0.021993 
det (3) = 0.164545 
det (3') = 2.00091 

A crude physical conclusion might be that, if the smallest effective action 
(— logdet/2) is preferred, on the basis of some unspecified dynamics, then the 3' 
twisting stands out as the lowest, for this artificial situation. The symmetry is 
broken from U(3) to U(l). 

5 The irreps are labeled by their dimensions. The subscript s stands for 'spinor'. 
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6. Application to thermodynamics on factored spaces 



In [21] a related analysis was made of Witten's symmetry breaking mechanism 
induced, as a toy model, by fluxes in geometries involving factors of the three- 
sphere. The Casimir energies were evaluated in the different gauge vacua and closed 
forms were found for lens spaces, the fields being in the adjoint representation of 
Q, = SU(n). The space-time was a factored Einstein universe, TxS 3 /r'. In [26,27], 
the system was put at a finite temperature, and symmetry restoration for increasing 
temperature analysed. The object of present interest, the effective action on S 3 /r', 
enters into this scheme through the high temperature expansions of thermodynamic 
quantities such as the free energy. 

General asymptotic series are given in [28] in terms of the coefficients of the 
short-time expansion of the heat-kernel, apart from one term which stems from 
the zero mode on the thermal circle. To a constant factor, this term is Z'(0, r", a). 
Furthermore, for the case of conformal propagation in four-dimensions, and scalar 
fields, the heat-kernel expansion terminates with the first volume, or Weyl, term, 
on fixed-point-free factors of the three-sphere, as already stated. The high temper- 
ature expansions, [28], of the, here finite, free energy, internal energy and entropy 
then simplify drastically to, (for a real scalar), 

E ~ r ' (10) 
s~^r 3 + ^'(o,r',o), 

respectively, with the error being exponentially small. 

These expansions have been rederived recently by Asorey et al, [29]. This 
reference also contains expressions for the complete thermodynamic quantities. 

Thermal theory on factors of the three-sphere had been set up by Kennedy, 
passim, [30], and was pursued in more detail by Unwin, [31]. 

As noted in [32] and [1], knowledge of the degeneracy generating functions 
allows one to write down the free energy, say, with minimum fuss, in terms of the 
cylinder kernel, K 1 / 2 , as, 



1 00 1 



m 

771=1 
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where E is the zero temperature internal energy or vacuum energy, or Casimir 
energy. 

In particular for twisted g-lens spaces, using (6), 



, , _ , , 1 1 d cosh(qm/35) 

m=l 



which is a way of rewriting the standard statistical sum over states. 

The twisted Casimir energies are effectively available in [1]. Equation (6) for 
the twisted cylinder kernel should be compared with equation (46) in [1]. The 
integer 'degrees' there, here transcribe to the reals, 

5 = q5, 6 1 =q/2, S 2 = l, 

and the Casimir energies are given in [1] equn.(52) 6 . Thus, 

(7 - 1205 2 + 2405 4 ) + 40(1 - 12S 2 ) q 2 + 112 



2880 q 



12) 



I have multiplied by two in (11) and (12) to allow for a complex field £ U(±)- 

Everything in (11) is explicit and computable. In figure 4, I plot the free 
energy versus the temperature, T = 1//3, for an untwisted field for q = 1,2, 3 and 
in figure 5 the free energy on the 4-lens for several twistings, i.e. fluxes. The high 
temperature behaviour, (10), can just about be seen. 

Fig. 4. q-lens free energy 
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" The S in the numerator should read Sq. 
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Fig. 5. 4-lens twisted tree energy 
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7. Two— sided quotients 

The same basic procedure can be applied to the inhomogeneous quotients, in 
the simplest case to the general lens spaces, L(q; X 1 ,X 2 ). The only difference is 
that, when when the degeneracy is being found, the necessary group average has 
to be done element by element. Without going into the standard details of the 
construction of the lens space, the elements of r" are determined by the angles a/3i 
and /3 2 , 



where p, = 0, . . . , q — 1 , labels 7. v x and u 2 are integers coprime to q, with X t and 



By an appropriate selection of a q-th root of unity, it is possible to set v x = 1, 
i.e. X x = 1, without loss of generality. Any pair, (z/ l5 z^ 2 ), can be reduced to (1, v) 
by multiplying through by the mod q inverse of u 1 . 



The simple, one-sided lens space, L(q; 1,1), corresponds to setting v = 1 so 
that 6 L = 0,6 R = 2-Kpjq. 



Pi _ Wi 
2tv q 



P2 _ V v i 
2tv q 



(13) 



A 2 their mod q inverses. 7 



Some authors label the lens space, equivalently, as L{q\ v ly v> 2 ). 
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The degeneracy generating function is, [1], 8 
^ «St=Z cob ft- cos & 



t(i t2) ^£(i+ t2 - l2tc °^i)( 



+ t 2 - 2tcos/3 1 7 VI + t 2 - 2tcos/? 2 



(14) 



with (13). The group sum over p has to be left until last. 

Again putting t = e~ T , the 4d-conformal cylinder kernel this time is, 



K 1 / 2 ( \ x \ SmnT 1 / 1t;l 

(r; A 1} A 2 ) - ^ _ ^ ^ )(coghr _ cog ^ , (15) 

with (13). When K 1 / 2 (r) is antisymmetric in r, as here and earlier, 9 the familiar 
contour routine can be followed leading to the general expression for Z'(0, a), 

i r°° /k 1 / 2 (Y n ) \ 

Z'(0,a) = -J dxRe\^ cosharj , r = x + iA , (16) 

(of which (4) and (7) are examples). A lies between and the imaginary part of 
the first pole of the cylinder kernel off the real r axis. 

As explained above, (13), it is sufficient to set u l = 1 and u 2 = v with v coprime 
to q. Then we have the corresponding — log det , 

rvi , \ 1 f°° , t-. ( sinhr coshar \ 

Z ^ a ^=q2^J ^ Re U(coshr-cosA)(coshr-cos^)J' (17) 

with < A < 2ir/q. This is easily computed. 

It is convenient to make the cyclic order, q, an odd prime so that all v from 1 
to q — 1 are relevant. I refer to v as the 'twist' of the manifold (as opposed to 5, 
the twist of the bundle which, in this section, is untwisted). 

The logdet for q = 29 and 1 < v < 28 is plotted in figure 6 for the 4d-conformal 
Laplacian. The effects of the isomorphisms between lens spaces can be seen. This 
is the same case figured in [2] computed by a different method so a comparison can 
be made. The results for the 3d-conformal Laplacian are very similar. 



° Higher dimensional spaces can be treated similarly e.g. [33,4]. A few remarks are in section 7. 
9 The antisymmetry arises essentially from the harmonic factor, (1 — t 2 ), in the generating 
function. 
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Fig. 6. logdet on twisted 29-lens space 
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The global finite temperature considerations of section 6 can be directly ex- 
tended to the inhomogeneous quotients. The Casimir energies would again be re- 
quired to find the analogue of (11) for the total free energy. Analytical results for 
these are contained in [34] although, to be in keeping with the spirit of this paper, 
a simple numerical summation of the group averages would suffice. 



8. Higher dimensions 

The notion of lens space extends to higher odd-dimensional sphere quotients, 
S 2e_1 /Z g (e.g. [35,9]). If q = 2e — 1 is an odd prime, they are catalogued by e 
integers, v i {% = 1 to e), the twists, coprime to q. The cylinder kernel for the 
Laplacian, conformal in 2e dimensions, is then, 

,v Vi , ) ^gn I _ i _i^, (18) 

where the angles P i are given by, 

2ttv 



(3 

and the previous numerical procedure is effected without difficulty. If one just 
requires a number, this approach is easier than the more analytical routes leading 
to Hurwitz ^-functions, say. 
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9. Minimal coupling 



The value a = 1 gives minimal coupling. There is then a zero mode (/ = 1) 
which shows up as an infra-red divergence as r tends to infinity in the relevant 
integrals such as (4). Removing the zero mode is conventional and amounts to sub- 
tracting e~ T from the cylinder functions. This however destroys their antisymmetry 
and the contour method used earlier is impossible. Instead I employ a different ap- 
proach which goes back to the early days of asymptotic analysis, see [36]. 

To begin with, I explicitly remove the (single) zero mode from the generating 
functions (cylinder kernels) by defining the subtracted functions, 

K 1/2 (r) = K 1 / 2 (r)-e- T 

H(r) = H(r)-e- T (19) 

z( s , r', a) = z( s , r', a) - 1 , 

(1 — a z ) s 

and note that Z(s, F', 1) is the minimal coupling ^-function by definition. 

Differentiating the basic (^-function definition (1) twice with respect to a 2 yields 
a higher resolvent, 



(7 \ 2 oo 
' ;— 9 



B ' (r<) (20) 



The sum converges at s = and so 



( d V zVn r' ^ 



1=2 v ' 



a 2 ) 2 



The idea is to integrate this quantity twice, determining the constants of integration 
from the reference point a = since Z(s, T', 0) = Z(s, T', 0) — 1, is known, [l]. 10 
First one has, 

°° T-% roo / \ 3/2 

E(^=^/ « r ; wy i ^ {aT) - (2i> 

I now use the fact that, (cf (3)), 



10 Any reference point could be chosen. This seems a convenient one, numerically. 
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to perform an integration by parts (the endpoint contributions vanish) , which gives, 

3/2 



= v /?£Vff W (£) 3/2 / 1/2 (ar) 



3/2 

(22) 



1 f°° —, . sinh ar 
dr r ii (r) 



2 Jo a 
A first integration with respect to a 2 is now easy to perform, yielding 

7 />0O 7 

z'(0,r',a) = drid(r) (cosh ar - 1) + — ; ,l'(0,r» 
Jo da 2 



da 2 



(23) 



a=0 

The constant of intregration (the final term) is found from the standard rela- 
tion, 

d °° D (T'l 

^2 ^ F '> «) = s E ( /2 _^ 2) i+i = + 1. F '' «) . (24) 

and one obtains for (23), 
d 



da 2 



/•OO 

z'(0,r',a) = / dr H(r) (cosh ar - 1) + Z(l,r',0) 

/•OO 

= / dr77(r)(coshar-l) + Z(l,r / ,0)-l, 
Jo 



(25) 



using (19). I evaluate the constant later. 

A final integration with respect to a 2 is required, and produces, 

/ f°° dT 

z'(0,r',a)= / — 77(r)(2ar sinh ar - 2 cosh ar + 2- a 2 r 2 ) + 

Jo T ' 

a 2 (z(i,r',o)-i)+z'(o,r',o), 

again using (19). The last term can be found numerically from the method in this 
paper, e.g. (4), as can the penultimate term. To be in harmony with the numerical 
viewpoint, the contour technique of [13] allows the values of Z{n,Y',a) (n G Z + ) 
to be obtained easily. For example, I find, after partial integration, for < a < 1, 

/•oo 

Z(l,T',a)= dx Re (H(t) cosh ar) , r = x + iA . (26) 
Jo 

As before, it is sufficient to consider the cyclic quotient, V = Z g , for which 
H(t) is given in (3) and so (r = x + iA, < A < 4/r/g), 

r,,. \ 1 f°° , ^ coth(or/2) cosh ar . 

Z(l,q,a) = - dxRe ^-Hr . 27 

v ' y ' ' 2 L sinhr v ' 
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Finally from (25), the minimal case is got by setting a = 1 to give my compu- 
tational formula for the minimal — logdet on homogeneous lens spaces, 

»- f % - -i — » - ^ (28) 

+ Z(l,g,0)-1 + Z'(0,g,0). 

One sees that the method of back integration has produced a systematic method 
of subtraction that simultaneously tempers the UV and IR divergences. 

The value a = gives the 4d-conformal operator and, although not needed, 
closed form expressions for the relevant quantities can be found in [1] . For example, 
for even-ordered lens spaces, one has the simple result, 

q-1 

Z(l, 2q, 0) = ^J2 cosec (^/S) • ( 29 ) 
q P =i 

This agrees numerically with (27) from which, of course, it can be derived. 

The case q = 2 provides a simple check. The quotient manifold is projec- 
tive three-space, P 3 . Periodicity restricts / to being odd, [17], and the ^-function 
rewrites to the difference of those on the three-sphere for minimal, W, (a = 1) and 
3d-conformal, 'c', (a = 1/2) propagation. Precisely, 

c^(«) = c53^)-2 2 - 2s a3(«), 

so that, 

C?S.'(0) = CF"(0)- 4^.(0), (30) 

where I have used the fact that C53 (0) vanishes, reflecting the closed nature of S 3 
and the absence of a zero mode. 

Both values on the right-hand side of (30) for the full sphere were evaluated an- 
alytically a long time ago giving Cp3 (0) = —0.695171, agreeing with the quadrature, 
(28). 11 

There do not seem to be many published numerical values for lens space deter- 
minants. Nash and O'Connor, [3], give some quite involved expressions with their 
corresponding plots, which seem to agree with my results. Figure 7 below gives an 
indication of my values, which, as an aside, suggest an inflection. For comparison, 
the upper curve is for the 4d conformal Laplacian (see fig.l), i.e. the last term in 
(28). 



11 Simplifications arise for this case. For example Z(l, 2, 0) vanishes. 
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Fig. 7. Minimal logdet on homogeneous q-lens 
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Using a = as a reference point is related to the method of binomially ex- 
panding the ^-function, (1), in powers of a 2 , utilised in [37]; see [38]. 

10. Discussion 

Relatively simple expressions allow the determinants on any homogeneous quo- 
tient of the three-sphere to be computed for any mass and for any twisting by an 
internal symmetry group. The inhomogeneous case has been treated in less detail. 

I have given some examples and plotted a few graphs, the physical significance 
of which is somewhat moot, and are meant only to illustrate possibilities. It is 
hoped that the numerical approach will be useful in confirming analytical results. 

I have discussed only scalar fields, but the extensions to spins half and one are 
straightforward, cf [1]. 

A different technique has been used for minimal coupling which involves back 
integrating the expression for a higher resolvent. This method could also be used 
for the other couplings although then the reference a = values would need to be 
found independently. 

Another set of quotients that lend themselves readily to similar treatment are 
the flat ones, classified by Hantsche and Wendt, on which a certain amount of work 
has been done. I draw attention now only to the early analysis of Unwin, [31,39], 
on the finite temperature field theory and resulting thermodynamics. Earlier, zero 
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temperature calculations can be found in [17,40] and de Witt, Hart and Isham, [41]. 
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